Boson pairing and unusual criticality in a generalized XY model 
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We discuss the unusual critical behavior of a generalized XY model containing both 27r-periodic 
and 7r-periodic couplings between sites, allowing for ordinary vortices and half- vortices. The phase 
diagram of this system includes both single-particle condensate and pair-condensate phases. Using 
a field theoretic formulation and worm algorithm Monte Carlo simulations, we show that in two 
dimensions it is possible for the system to pass directly from the disordered (high temperature) 
phase to the single particle (quasi)-condensate via an Ising transition, a situation reminiscent of the 
'deconfined criticality' scenario. 



More than 25 years ago, Korshunov pQ and Lee and 
Grinstein [2| discussed the statistical mechanics of cer- 
tain generalizations of the familiar XY model. XY mod- 
els are the simplest systems capturing the rich physics of 
vortices, and these topological defects are well known to 
govern the critical behavior in two dimensions [3 . The 
generalizations discussed by the above authors give rise to 
half vortices, about which the XY order parameter winds 
by 7r, connected by strings with a finite tension. This in 
turn leads to a far richer phase diagram (see Fig. [I]), 
whose essential elements were verified by numerical sim- 
ulation [4]. 

Over time there has been considerable interest in iden- 
tifying systems where the physics of half vortices and 
strings plays a role, ranging from nematic liquid crystals 
[21 [5] to the A-phase of 3 He [HEj, to spinor Bose conden- 
sates 7, 8 . Indeed, half vortices have now been directly 
observed in exciton-polariton condensates [SJ. In recent 
years a large amount of activity has focused on one par- 
ticular candidate: a gas of attractive bosons. It can have 
two distinct superfluid phases: an atomic superfluid of 
bosons and a molecular superfluid of boson pairs, with 
the latter supporting half vortices [TU1 E] • Such a sys- 
tem is in general likely to be unstable to collapse, but a 
possible resolution of this difficulty is to harness three- 
body loss to project out triple occupancy of each site of 
an optical lattice [12]. This idea led to a resurgence of 
interest in the problem [T3TU7| . 

We have looked anew at the phase diagram of this type 
of system. Previous studies found that one could cither 
pass directly from the normal state to the atomic conden- 
sate at low temperatures, or first into a molecular conden- 
sate with no single-particle long-range order, and thence 
to the atomic condensate in an Ising transition. Remark- 
ably, we find that in two dimensions there is a region of 
the finite temperature phase diagram where the normal- 
to-atomic superfluid transition is of the Ising, rather than 
the Kosterlitz-Thouless type (see Fig. [I]). Such uncon- 
ventional behavior is reminiscent of the 'deconfined crit- 
icality' scenario [18) . The two have a common origin 
in that the expected proliferation of point-like defects 
is suppressed by critical fluctuations. Our prediction is 
further supported by Monte Carlo simulations using the 
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FIG. 1. Schematic phase diagram of the model Eq. (2f. Solid 
lines are continuous transitions of the type indicated; the 
dashed line is first order. The dotted line indicates the region 
where the nature of the transition is unknown. The Renormal- 
ization Group analysis developed here applies in the vicinity 
of the solid dot. Inset: an illustration of a half vortex pair 
and the associated string. 



worm algorithm and a novel method based on conformal 
field theory to identify an Ising transition in a system 
with non-Ising degrees of freedom. We stress that our 
result applies to the finite temperature phase diagram 
of many of the above mentioned systems in 2D (see e.g. 
[S[ [SI EH E])) as well as to the zero temperature phase 
diagram of the ID system in Ref. [TT] . 

The simplest generalization of the XY model with the 
requisite physics has the form 



ff = -\J(l-A)cos( 

{ij) 



9j) + A cos (20* -20J-)] (1) 



where the angular variables B. L : < Bi < 2tt are defined 
on a square lattice and (ij) represents nearest neighbor 
pairs. In quantum mechanical language, the two terms 
correspond to the hopping of single bosons and boson 
pairs respectively. A = recovers the usual XY model, 
which allow for the existence of point-like vortices around 
which the phase increases by a multiple of 27r. A = 1 cor- 
responds to an XY model with half the periodicity, and 
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thus half vortices of 7r winding, but otherwise identical 
properties. For A > 1/5 the intersite energy develops a 
metastable minimum at \6i — 0j\ — 7T. Pairs of half vor- 
tices are connected by a 'string' where 6^ — 8j jumps by ir 
and for A < 1 this string has a finite tension (see Fig. [T]). 

When this tension is large (small A) half vortices are 
bound together irrespective of their topological charge, 
so that only integer vortices may exist freely. Thus here 
the model displays the familiar Kosterlitz-Thouless (KT) 
transition. When the tension is small (A close to 1), a 
KT transition of the half vortices can occur (|KT). After 
this transition half vortices are bound together in pairs, 
so that the strings connecting them form closed domain 
walls, which disappear at a lower temperature as the ten- 
sion overcomes their entropy. This second transition is 
of the Ising type. 

The resemblance of strings to Ising domain walls that 
can terminate on half vortices suggests that the operator 
inserting a half vortex includes an Ising disorder opera- 
tor [19] . In the high temperature phases (corresponding 
to the disordered and pair superfluid phases in Fig. [I]) 
the strings are not confining, equivalent to long range or- 
der for the disorder operator. The half vortices then can 
drive a KT transition by the familiar mechanism. How- 
ever, along the Ising critical line the expectation value 
of the disorder operator vanishes. The resulting critical 
fluctuations of the strings suppress the proliferation of 
half vortices, so that a direct Ising transition between 
the disordered and superfluid phases can occur. 

Villain model. This qualitative picture is borne out by 
our analysis of a particular microscopic model. Following 
the original references [TJ |5] , we study a variant of the 
model Eq. ([I]) with partition function 

2 = \{f II M*«6) + e- K w v (0 ab - 7T)) (2) 



Integrating out the angular variables leaves 
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with J* = J- 1 . K -t oo is the usual Villain model, with 
a KT transition at Jyuiain ~ 0.75 . . . [20j[2T], while K = 
corresponds to a n periodic Villain model, allowing free 
half vortices, with a |KT transition at J = 4Jyin a in ~ 3. 
Finite K gives the strings connecting half vortices a finite 
domain wall energy. When J oo the phase differences 
between neighboring sites are restricted to or ir and 
the usual square lattice Ising model with transition at 
Rising = ln(l + v2) ~ 0.881 ... is recovered. 

We map to a generalized height model using the sec- 
ond representation of the Villain potential in Eq. Q. 
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where and K satisfy the Kramers-Wannier duality 



relation sinh sinh K = 1. Here the variables 



live 



on the links of the lattice. Due to their vanishing lattice 
divergence, they may be thought of as currents describing 
the world lines of bosons. 
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where the variables \ii = (—l) hi . The discrete Gaussian 
model given by the first term of Eq. ([5| has a roughening 
transition from a smooth phase at small J (high temper- 
ature) to a rough phase at large J, corresponding to a 
superfluid state of the bosons. The new term in K t al- 
lows for the existence of a second rough phase where the 
heights are predominantly even or odd, with the currents 
riij even, corresponding to a pair superfluid. 

The Ising and Gaussian parts are disentangled by writ- 
ing hi = 2hi + [fj,i — 1)/2 and identifying the fii as the dis- 
order variables of an Ising model, dual to the usual spin 
variables. We move to continuous hi by introducing the 

J-functions IL E P <^ ~ p) = Ui E g (?™ qh ■ An effec- 
tive sine-Gordon description is arrived at in the standard 
manner |22j by introducing an effective coupling (vortex 
fugacity) z p for each harmonic. After a shift in the inte- 
gration variables hi — > hi — (/ij + l)/4 and retaining only 
the q — ±1, ±2 terms we arrive at 

z = E /n^ ex p(E^+E^) . 



Vi = z\ /ii cos(27r/ij) + zi cos(47r/ij) . 

The nearest-neighbor interactions Sij are those of decou- 
pled discrete Gaussian and Ising models. The two terms 
in the potential Vi describe respectively the half and inte- 
ger vortices, making precise how the half vortices couple 
the Ising and Gaussian degrees of freedom. 

Renormalization Group (RG) analysis. The limits 
shown in Fig.[l]provide the skeleton of the phase diagram 
for this model, but to understand what happens when 
the transitions approach each other requires an analy- 
sis of the coupling between the Gaussian and the Ising 
degrees of freedom. Since we are concerned principally 
with the nature of the transition in this region, this is 
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most conveniently accomplished in the field theory limit, 
where we have a sine-Gordon model and an Ising model 
coupled by the perturbation 

H i_ v = z\ J dx/i(x) cos[27r/i(x)] . (6) 

In the high temperature (low K) phase of the Ising 
model, the disorder operator /i(r) acquires an expecta- 
tion value. Here Hi_ v becomes equivalent to the cosine 
potential of the sine-Gordon model. When z\ = Z2 = 0, 
the operator cos(2irnh) has dimension 7rn 2 /(4J„), and 
so is relevant for J < 8/ (7m 2 ). Thus in the Ising high- 
temperature phase, the half-vortices and vortices drive 
KT transitions at J = 8/tt and J = 2/ir respectively. 
(These should be understood as the renormalized values 
of J: the transitions occurs at a larger value of bare J.) 
This accounts for the fourfold increase in the jump of the 
superfluid density at the |KT transition relative to the 
usual one . 

The unusual behavior occurs along the critical line of 
the Ising model. Here /i(x) no longer has an expectation 
value, but only critical fluctuations obeying (/i(x)/z(y)) ~ 
|x — y I -1 / 4 . The scaling dimension of /i(x) cos[27r/i(x)] is 
I + -rr > which does not become relevant until J = ^ , a 
smaller value than for proliferation of half vortices in the 
disordered phase. This suggests the scenario depicted in 
Fig. [T] the Ising transition persists after it has met the 
^KT transition. 

The RG equations give more insight into these transi- 
tions. Defining the deviation from Ising criticality to be 
k = K — K c , at Z2 = we find to second order in k and 
zi [H] 
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where I = log(£/£o), with £ being the coarse-grained 
length scale. In Fig. [2] we show two sections of the flow in 
the k — zi plane, for values of J above and below J c = ip. 
When J > J c , the Ising fixed point at the origin is stable 
to the Hi_y perturbation. Note, however, that this is a 
dangerously irrelevant perturbation, with z\ growing at 
negative k (higher temperature) reflecting the prolifera- 
tion of the half vortices. In contrast, when J < J c , the 
Ising fixed point is unstable. The apparent fixed point 
at finite z\ is in fact a separatrix along which J flows to 
zero: crossing this separatrix presumably corresponds to 
a first order transition. 

The RG equations Eqs. ^ give are valid in the vicinity 
of the Ising fixed point. In particular, our conclusion 
regarding the dangerous irrelevance of Hi_ v applies close 
to the critical value J c = k& where this perturbation 
becomes marginal. When J > — we know that Hi v is 

o 7T 2 -V 

irrelevant even deep in the disordered phase. 
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FIG. 2. (Left) RG flow in the k - Zi plane for J = §i. 
Here the fj,cos(2irh) perturbation is irrelevant, even though 
(fj) cos(2-7rft) is relevant. Note that the Ising fixed point at 
k = 0, Z! = is stable. (Right) RG flow for J = ff . The 
Ising fixed point is now unstable. 
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FIG. 3. Phase diagram near the meeting of the Ising (O) 
and -^KT (x) transitions determined by numerical simula- 
tion. The transitions are determined as described in the text 
with the Ising line stopping at the last point where we are 
confident of the nature of the transition. The dotted line is a 
sketch of the (presumably) first order line connecting to the 
usual KT transition (see Fig. [T|. 

Numerical simulations. We have tested the above cal- 
culations using Monte Carlo simulations of the model 
Eq. ([2| based on the worm algorithm [25] . in the formu- 
lation given in Ref. |26) . Both single and double worms 
are needed to accurately simulate the paired phase. The 
worm algorithm directly simulates the partition sum in 
term of currents, Eq. Q. The phase diagram resulting 
from our simulations is illustrated in Fig. [3j clearly show- 
ing the persistence of the Ising transition after the Ising 
line meets the |KT line. 

A particular advantage of the worm algorithm is that 
it simultaneously simulates all sectors of the model de- 
fined with periodic boundary conditions (i.e. on a torus) . 
A sector is specified by a pair of integers W = (W X7 W y ) 
giving the winding of the current loops around each circle 
of the torus; in terms of heights, h(x+L,y) = h(x,y)+W y 
and h(x,y + L) = h(x,y) + W x . To locate the KT and 
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FIG. 4. The ratio Eq. (JsJ) of partition functions is used to 
locate the Ising transition from the crossing point for different 
system sizes. Here O, +, □, x symbols are for systems sizes 
L — 20, 40, 80, 160 respectively, J = 2.8, and the (small) error 
bars are included. The dotted line corresponds to the pure 
Ising value fc = .3729. 



|KT phase boundaries, we follow the method of Ref. [57] 
and exploit the fact [28] that the superfluid density (he- 
licity modulus) is T = |^(W 2 ) (we restore temperature, 
which provides the energy scale). 

To locate the Ising transition, we utilize the partition 
functions in the various sectors. The non-universal bulk 
contribution to the free energy is independent of sector, 
and so the ratio of partition functions in different sectors 
should be a universal property of the critical point. Thus 
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becomes independent of system size, and so the curves 
for different sizes plotted in Fig. [3] cross at the transition. 
We extract the correlation length critical exponent by as- 
suming the form £(L, J, K) = /(L/£) with the correlation 
length £ ~ kT v applying near the critical point. The best 
scaling collapse (Fig. [5]) is achieved with v = 1.00 ±0.02, 
consistent with v — 1 for the Ising model. 

In fact, the critical ratios themselves are known exactly 
from conformal field theory [33] . Even though the Gaus- 
sian and Ising actions decouple at z\ = 0, the fields still 
satisfy a nontrivial 'gluing condition'. If W x (say) is odd, 
then the /x^ satisfy antiperiodic boundary conditions in 
the y direction, and h(x, y + L) — h(x, y) is half integer. 
In the field theory limit, we have at the decoupling point 
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where 'P' ('A') denotes the (anti-)periodic boundary 
conditions in the Ising part, and (1/2) denotes 
(half) integer winding in the bosonic part. For the 
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2tf 2 (0,e- 7r ) / %(0,e- 7r ) = 2.68179... , where the &j(z,q) 
are the standard Jacobi theta functions. The Gaussian 
contribution to C is tf 3 (-7r/2, e" 7r2j / 2 )/i9 3 (0, e^ 2 ' 7 / 2 ), 




FIG. 5. C plotted vs. k,L 1/v . For the same L as in Fig. [ZJ the 
best collapse is achieved for v = 1.00 ± 0.02 



equal to unity to the fifth significant digit at J = J c , 
and increasing to one in the region of interest J < J c . 
Thus the ratio £ at the Ising critical point is close to 
C c = .3729. In Fig. [4] the values of £ at the crossings are 
indeed approaching £ c with increasing system size. 

In conclusion, we have shown that a two dimensional 
XY system supporting half vortices and strings has un- 
usual critical behavior driven by the interplay of these 
two types of defects. The same physics is expected to 
play a role in the (1 -I- 1) -dimensional quantum problem 
(see Refs. [32] [17]), where there is the additional freedom 
of lattice filling to explore. In three dimensions (or the 
2 4- 1-dimcnsional quantum problem) the analog of the 
coupling between the half vortices and the disorder vari- 
ables of the Ising model is an Ising gauge charge carried 
by half vortex lines. The consequences of this coupling 
will be explored in future work. 
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